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For many years a spin-state transition at T ~ WOK and insulator - metal transition 
(IMT) at TiMT ~ QQQK in LaCoOa remains a mystery. Small low-spin - high-spin spin gap 
As = E (HS) — E (LS) ~ WOK results in the spin-state transition. The large charge gap 
2Ea ~ 2300A' {Ea is the activation energy) vs. As and Timt implies that LaCoOa is not a 
^~^. simple narrow-gap semiconductor. Here we explain both the spin-state and IMT on the same foot- 

ed ' ing. We obtain strong temperature dependent band structure in LaCoOs by the LDA-I-GTB method 

that incorporates strong electron correlations, covalence and spin-orbital interaction exactly inside 
the CoOe cluster and the interclaster hopping between different multielectron configurations by 
perturbation theory for Hubbard X-operators. 



PACS numbers: 72.80.Ga, 71.30.+h, 71.70.Ch, 75.30.Wx 



The perovskite-oxide LaCoOa has been studied intensely for many years due its unique magnetic properties and 
related insulator- metal transition (IMT) [l|, I4I . A gradual appearance of the paramagnetism above SOX from the 
nonmagnetic ground state in called the spin-state transition. Goodenough was the first who suggested that in- 
j^ , stead of the Hund's rule dictated high spin (HS) S = 2 the strong crystalline electric field results in the low spin 
(LS) S — Q state for cf' configuration of the Co'^^ ion, and the energy difference is very small with the spin gap 
T • A5 = E {HS) — E [LS) ~ lOOX. The thermal population of the HS state provides the sharp increase of the magnetic 
■ O ' susceptibility x with a maximum around lOOiC. The nature of the excited spin state of Co'^"'" above the singlet ^^lo 
5 . has been under debate (see for a recent review [3[). Besides the original ^T2g HS state with the ifgCg configurations [3| 
fj ■ there were many indications on the intermediate spin (IS) S = 1 ^Tig state. The two stage model has been proposed 
with the LS-IS transition at \QQK and IS-HS transition at 550 — GOOivT IJ,|5J. Recent electronic spin resonance (ESR) 
[6|, X-ray absorption spectroscopy (XAS) and X-ray magnetic circular dichroism (XMCD) [7| experiments prove that 
";> ' the lowest excited state is really the HS. Nevertheless the ^T2g term is splitted by the spin-orbital interaction in 
CN| \ the low energy triplet with effective moment J = 1, and higher energy sublevels with J = 2 and J = 3 [8!]. The 
C^ I large difference between the spin excitation gap A5 and the charge gap given by the activation energy for electrical 
l^ ' conductivity Ea ~ 0, leV at low T indicates that LaCoOs is not a simple band insulator [9[. The second shallow 
maximum in x near 500 — 600K is often related to the insulator-metal transition. Surprisingly for the IMT electrical 
1^ conductivity a does not seem to show any noticeable anomaly at this temperature [9|. Moreover the discrepancy 
between the large charge gap 2Ea ~ 2300K and the insulator-metal transition t emp erature Timt ~ GOOif implies 
that the IMT cannot be simply argued in terms of a narrow-gap semiconductor [lO|. Here, we solved this problem 
by calculating the electronic band structure in the regime of strong electron correlations. We consider electron as the 
linear combination of quasiparticles (QP) given by excitations between the different multielectron configurations ob- 
tained by exact diagonalization of the CoOe cluster. With the Hubbard operators constructed within the exact cluster 



5—1 ' eigenstates we can calculate the QP band structure for the infinite lattice. The QP spectral weight is determined by 
. . . the occupation numbers of the local multielectron configurations. We find that the thermal population of different 
sublevels of the ^T2g HS term splitted by the spin-orbital interaction results both in the spinstate transition and also 
in some new QP excitations. Of particular importance is the hole creation QP from the initial d^ HS into the d^ HS 
term, this QP appears to form the in-gap state inside the large charge-transfer gap Eg « l,5eV. The intercluster 
hopping transforms this local QP into the in-gap band that lies just under the bottom of empty conductivity band 
and provides the insulating gap 2Ea ~ 0, 2eV at T = lOOK. It bandwidth increases with T, and verlaps with the 
conductivity band at T = Tjmt = 587K resulting in the IMT. Hence our approach allows to treat both the low T 
spin-state transition and the high T IMT on the same footing. 

LaCoOa as well as other strongly correlated oxides is a difficult problem for the ab initio band theory. The LDA 
calculations fllj incorrectly predict a metal for paramagnetic LaCoOa. Various methods have been applied to study 
effect of correlations on the LaCoOa electronic structure: LDA-I-U or GGA-I-U [12l4l5l |. dynamical mean-field theory 
[16| . Recent variational cluster approximation (VGA) calculation [l7| based on the exact diagonalization of the 



CoOg cluster gives a reasonably accurate description of the low temperature properties: the insulating nature of the 
material, the photoelectron spectra, the LS-HS spin-state transition. The main deficiency of the VGA is the failure 
to reproduce the high temperature anomalies in the magnetic and electronic properties associated with the IMT. The 
exact diagonalization of the multielectron Hamiltonian for a finite cluster provides a reliable general overview of the 
electronic structure of the correlated materials [18|. To incorporate the lattice effect several versions of the cluster 
perturbation theory are known [l9|, |20|. To calculate the band dispersion in the strongly correlated material one 
has to go beyond the local multielectron language. The natural tool to solve this problem is given by the Hubbard 
X-operators X^'^ = \p) {q\ constructed with the CoOe cluster eigenvectors \p) at site Rj. All effects of the strong 
Coulomb interaction, spin-orbit coupling, covalence and the crystal field inside the CoOe cluster are included in the 
set of the local eigenstates Ep. Here p denotes the following quantum numbers: the number of electrons (both 3d 
Co and p of O), spin S and pseudoorbital moment I (or the total pseudomoment J due to spin-orbit coupling), the 
irreducible representation in the crystal field. A relevant number of electrons is determined from the electroneutrality, 
for stoichiometric LaCoOs n = 6. In the pure ionic model the corresponding energy level scheme for dP Co^~^ has 
been obtained in [8|. Due to the covalence there is admixture of the ligand hole configurations d"+^L and d"^^L^ 



that is very well known in the X-ray spectroscopy 21| . Contrary to spectroscopy the electronic structure calculations 
require the electron addition and removal excitations. For LaCoOa it means the dP and d"^ configurations. The total 
low energy Hilbert space is shown in the fig. 1. Here the energy level notations are the same as in the ionic model [8| 
but all eigenstate contains the oxygen hole admixture due to the covalence effect. The calculation of the n = 5, 6, 7 
eigenvectors for CoOe cluster with the spin-orbit coupling and the Coulomb interaction has been done in 



22 



Following the generalized tight-binding (GTB) method ^] and its ab initio version LDA-I-GTB [2j| we consider the 
electron annihilation operator at site /, orbital state A (Co d or O p) and spin projection a as the linear combination 
of the Hubbard operators 



afXa 



= J2^,^in)Xf. (1) 



Here n numerates different pairs of indexes in XP'' = \p) {q\ to describe the excitation from the initial state \q) to 
the final one \p). The matrix element 7act("-) is calculated straightforwardly as the eigenvectors {\p)} are known from 
the exact diagonalization of CoOg cluster. In the X-operator representation the initial multielectron and multiorbital 
Hamiltonian has the same operator structure as the Hubbard model 

fp fgnm 

Here Ep is the local eigenstate. All intercluster Coulomb interactions and hopping are included in the tfg matrix and 
are treated perturbatively. The electron Green function 



GAA'.(fc,c^) = {{akxa I a+y,))^ = ^ 7A.(n)7A'.(m)D™"(fc,c^), Df^^ = (^(^ 
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Its local part is found exactly 



GTAk,u)^Y.\^Un)?^. (4) 

n " 

Here r2„ ^ Ep{N +1)- Eg (N) is the QP energy, F{n) = (^^^) + (xj"^) is the fiUing factor. This factor provides 
zero spectral weight for the excitations between two unoccupied states. The non-zero contribution to the Green 
function requires a participation of at least one occupied state in the QP n = (p, q) . 

The electron removal spectrum determines the top of the valence band, the corresponding QP are shown in the 
fig. 1 by thin solid lines as the excitation from the ^Aig d^ singlet in the ^T2g d^ states with J = 1/2 and J = 3/2. 
There energies are 

nvi = E (d^ ^Ai) - E (d^ 2t2, j = 1/2) ,nv2^E (d^ ^ai) - e (d^ ^ts, j = 3/2) . (5) 

The bottom of empty conductivity band has the energy 

nc = E{d\^E)^E{d\^Ai). (6) 



All these bands have nonzero QP spectral weight. The intercluster hopping results in the dispersion, r2„ — ?► iln{k). 
In a simplest Hubbard- 1 approximation for the tfg Green function is given by 

D-\k,u:)^D,\oj)~i{k), (7) 

the band structure is obtained from the dispersion equation 

det \\Smn [E - nn)/F{n) - i'""(fc) || = 0. (8) 

In the LDA+GTB version the single-site electron energies as well as the interatomic hopping integrals of the p-d model 
has been obtained by Wannier function projection technique [25| . The LDA band structurejaf LaCoOa was obtained 
within linearized muffin-tin orbital (LMTO) basis set [26| using structural data for T — 5K [2^. The calculated bands 



are in good agreement with the previous results 11|. The bands resulted from the projection procedure for five Co 
d-orbitals and three O p-orbitals completely reproduce the initial LDA bands. The local coordinate system with the 
axis directing towards oxygen atoms forming the octahedron around Co atom was used. The intra-atomic Hubbard 
and Hund parameters U and J have been calculated within the constrained LDA method [28|. The QP LDA+GTB 
band structure corresponds to the charge-transfer insulator [29| with the gap Eg w 1.5eV (fig. 2) at T = OK. This 
gap value is rather close to the VGA gap [17| and the experimental value Eg « leV [lOi] . 



At finite temperature the thermal excitation over the spin-gap Ag into the J = 1 and over the gap A5 -f 2A into 
the J = 2 sublevels of the HS ^T2g state occurs. We take Ag — lAOK and A — 185K following6. Partial occupation of 
the excited HS states results in the drastically change of the QP spectrum. For T = OK excitations from the ^Ai dP 
singlet in the lowest ^Ai dP term were forbidden due to spin conservation (the corresponding matrix element 7„ = 0, 

see (IlJ), and the excitation from d^, J = l\ in Itf', ^^1) has nonzero matrix element (shown by dashed line fly^ in 

the fig. 1) but zero filling factor as the excitation between two empty states. For T ^ the filling factor for the fty-^ 



and r2y2 QP is non zero and is equal to the occupation number ni and ^2 of the states 
correspondingly. The energies of these QP are 



d^ J =1) and 



(f, J = 2 



^vi — E 



(d^ 'T2g, J=l)-E {d\ ^Ai) ,n*y2^E (d^, 5T2,,, J^2^-E {d\ "^ A^) . (9) 



The energies of these QP appear to be slightly below the bottom of the conductivity band, see DOS at finite temper- 
ature in the fig. 2. Thus we have obtained that these temperature-induced QP states lies inside the charge-transfer 
gap, they are the in-gap states. Similar in-gap states are known to result from doping in the high temperature super- 
conductors. The LaGoOa is unique because the in-gap states are induced by heating. The chemical potential lies in 
the narrow gap 2Ea ~ 0.2ey at T = lOOK between the in-gap states and conductivity band. 

From the GTB dispersion equation (|S]) it is clear that the in-gap bandwidth is proportional to the occupation 
numbers rii and 712 of the excited HS states. With further temperature increase the in-gap bands ^yi and £7^2 
become wider and finally overlap with the conductivity band ^c (fig- 2) at T = Timt — 587K. It should be clarified 
that the IMT in LaCoOa is not the thermodynamic phase transition, there is no any order parameter associated with 
the gap contrary to the classical IMT in VO2, NiS etc. 

Assuming that the carrier mobility is a smooth function of temperature we have estimated the temperature depen- 
dent conductivity as 

a(T) = f7oexp(-^,(T)/fcT) (10) 

with the temperature dependent Ea{T) given by the DOS calculations. The fitting parameter <tq was taken as 
Co = CTcxp {800K). The overall agreement of the calculated and experimental data from [9| is clear (see fig. 3). The 
deviations at the 300 — 2>^0K we believe stem from some other interactions omitted in our calculations. There is 
similar anomaly in the thermal expansion coefficient in the same temperature range 3(| , so may be the spin-phonon 
and electron-phonon interactions are responsible for this discrepancy. 

We have also calculated the average moment J = L + S* as Jav ~ \J'^) as function of temperature (fig. 4). Here 



^') = EM'^'h)(^"")- (11) 

n 

The expected for HS value Jav « 2 is reached only at T w lOQOisr. In the region of the spin-state transition at 
T « lOOK the value of Jav is close to 1. May be it is the reason why so many experimentalists have obtained the IS 



5 = 1 fitting their data by the state with definite spin. Temperature dependence of the magnetic susceptibihty has 
the maximum at T w lOOK similar to many previous works with LS-HS scenario |9l |31|. 

Thus, we find that a correct definition of the electron in strongly correlated system directly results in the in-gap 
states during the spin-state transition due to the thermal population of the excited HS states. Close to the spin-state 
temperature region the in-gap states determine the value of the activation energy Ea ~ O.lel^. Further temperature 
increase results in larger in-gap bandwidth and smaller Ea, and finally Ea ^ at Tjmt = 587 K. As concerns the 
weak maximum in the x(^) close to the IMT it may be a small Pauli-type contribution from the itinerant carriers 
above Tjmt- We emphasize that instead of rather large difference in temperatures of the spin-state transition 
(^ lOOK) and the IMT {600K) the underlying mechanism is the same and is induced by the thermal population of 
the excited HS states. 
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FIG. 1: The low energy part of the CoOs cluster Hilbert state for the electron numbers A'^e = 5, 6, 7. Terms with given Ne 
are the mixtures of d^', d'^^'^^L^ and d^^'^^L^ configurations. At T = only the N^ = 6 low spin term ^Ai is occupied, the 
Fermi- type excitations from this term that form the top of the valence band (d® — > d^) and the bottom of the conductivity 
band (d^ — >■ d^) are shown by the solid lines with arrow. The dashed lines denote the in-gap excitations with the spectral 
weight increasing with temperature due to the population the HS excited d® terms. 




FIG. 2: Single particle density of states at different temperatures. At T = OK LaCoOa is the charge transfer insulator with the 
gap Eg ~ 1.5el/. At finite temperatures the in-gap band appears below the conductivity band with the temperature dependent 
activation energy. At T = lOOA' Ea ~ O.Iel/. At T = Tjmt = 5S7K Ea = OeV, and above Timt the band structure is of the 
metal type. 




FIG. 3: The temperature dependence of the resistivity, solid line shows the experimental data [3| and dashed line is calculated 
from Eq.(4). 
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FIG. 4: The temperature dependence of the average local moment shows the spin-state transition from the nonmagnetic low 
spin state at T = OK to the paramagnetic high spin state at finite temperature. Only at T ~ lOOOii' the moment is close to 2. 
In the region of the spin-state transition at T ~ lOOK Jav ~ 1. 



